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Abstract 

We study the motion of a heavy tracer particle weakly coupled to a dense interacting Bose 
gas exhibiting Bose-Einstein condensation. In the so-called mean-field limit, the dynamics 
of this system approaches one determined by nonlinear Hamiltonian evolution equations. 
We derive the effective dynamics of the tracer particle, which is described by a non-linear 
integro-differential equation with memory, and prove that if the initial speed of the tracer 
particle is below the speed of sound in the Bose gas the motion of the particle approaches 
an inertial motion at constant velocity at large times. 
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1 Background from Physics and Equations of Motion 

In this paper we study the motion of a very heavy tracer particle coupled to a very dense, very 
weakly interacting Bose gas at zero temperature, which exhibits Bose-Einstein condensation. 
In an interacting Bose gas at positive density and zero temperature, the speed of sound is 
strictly positive. If the initial speed of the tracer particle is well below the speed of sound in 
the gas one expects that the motion of the particle will approach a uniform ballistic (inertial) 
motion at large times. A result in this direction has recently been established in a certain 
limiting regime (the "mean-field-Bogolubov limit") of the Bose gas in p]. In the present 
paper, we prove some results complementary to those in [2J for the same model. (For earlier 
results on related models, see also [3 E].) 

We consider a tracer particle of mass AM and a Bose gas of density Apo with two-body 
interactions of strength A~^k. In the so-called mean-field limit of the quantum system, where 
A — )■ cxD, see [H [I], we obtain the following classical equations of motion 



In Eqs. (11. ip and (II. 2p . Xf G M.^ and Pj G M.^ are the position and momentum of the 
tracer particle at time t, respectively, and at{x) is the Ginzburg-Landau order-parameter 
field describing the state of the Bose gas at time t, in the mean-field limit; W and are 
two-body potentials of short range, with of positive type, and g and n > are coupling 
constants. The interpretation of |af(a;)p is that of the density of bosonic atoms at the point 
X of physical space M^, at time t. The global phase of ctt is not an observable quantity. The 
symbol * in (II. 2p denotes convolution. 

Eqs. (II. ip and (II. 2p turn out to be the Hamiltonian equations of motion corresponding 
to the following Hamilton functional 



^t=^, Pt = -VxV{Xt)+g j dxV,W{Xt~x){\at{x)\^~^}, 
(x) = f-^A + gW{Xt - x)\ at{x) + kU* {la^l^ - ^j") (x) at{x), 



(1.1) 



(1.2) 



H{X, P, a, a) 



p2 




2M 
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The Poisson brackets are given by 

{X\X^} = {P„P,} = 0,{X\P,} = S], (1.4) 

and 

{J{x), J{y)} = 0, {«(x), a{y)} = -tS{x - y). (1.5) 

We impose the conditions that Vat is square-integrable in x and that \at\^ — ^ is in- 

tegrable. In the present paper, we also require that at{x) — )■ as |x| — )■ oo. This last 

condition explicitly breaks the gauge invariance, a^(x) — )■ e^*^a''(a:), where 9 is an arbitrary 
angle (independent of x and t). Given our boundary conditions on o;" at oo, it is natural to 
define a new function /3 by 

atix)=: .m+/3tix), (1.6) 



with f3t{x) — > 0, as |x| — )■ oo. The equations of motion then read 

Xt Pt = -VxV{Xt) +gj V,W{Xt - x) (^\/3t{x)\^ + 2^RePt{x)^ dx, 



;i.7) 



^^{x) = (^-^A + gWiXt - x)^ A(x) + ^W{Xt - x) 




The Hamilton functional giving rise to these equations is obtained from (11. 3p after inser- 
tion of Eq. (11. 6p . It is easy to see that, under rather weak assumptions on the potentials 
W and 0, Eqs. (11. 7p and (II. 8p have static solutions, and that if the external force vanishes 
[V = 0) they have "traveling wave solutions" , provided the speed of the particle is smaller 
than or equal to the speed of sound in the Bose gas; see |5l [2]. These solutions correspond 
to an inertial motion of the tracer particle at constant velocity, and the particle is accompa- 
nied by a "splash" (a coherent cloud) of atoms from the Bose gas. (Quantum mechanically, 
this splash corresponds to a coherent state of gas atoms and causes decoherence in particle- 
position space, which allows for an essentially "classical" detection of the particle trajectory.) 
If, initially, the speed of the tracer particle is larger than the speed of sound it emits sound 
waves - Cherenkov radiation - into the condensate, which causes friction. As a consequence, 
the particle loses kinetic energy until its speed has dropped to the speed of sound in the Bose 
gas. This phenomenon has been analyzed for a simple model (the B-model defined below) 
in [3]; (work on Cherenkov radiation in more complicated models is presently carried out). 

The following models are of interest (see [3, S]): 
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B -Model: k = (ideal Bose gas), — ?■ 0, see [1]. 
C -Model: K = 0,hnt 0; see [3]. 

E -Model: 2Kpo/g'^ := A =const., with g,K ^ ("Bogolubov limit"). 
G -Model: k > and ^ ^ 0. 

In this paper, we focus our attention on the E-model, with V = and 4>{x) = 6{x). The 
equations of motion then take the form 

Xt =^Pu Pt = v/p^i?e(V,iy^% A) (1.9) 

tk = - T^AA + XRef3t + v^W^^', (1.10) 
2m 



W^{x) ■.= W{X -x). (1.11) 



where 



The speed of sound is given by 

The Hamilton functional of the E-model is 



i/(X,P; /3,/3) := + — / \V[5\'dx + \l \Rel3\''dx + 2^ I WRe/Sdx. (1.12) 
2M 2m 



In the subsonic regime the equations of motion have a family of solutions describing 
inertial particle motion ( "traveling- wave solutions"). The existence and stability of these 
solutions has been studied in [2]. In this paper we extend the results of [2]. We consider 
initial conditions at time t = with the properties that /So is "small" in a suitable sense, 
i.e., the state of the Bose gas is close (or equal) to the ground state [3q = 0, and that the 

initial speed of tracer particle is well below the speed of sound -y/^. We prove that if the 

interaction potential W is sufficiently weak the particle will approach an inertial motion at 
a constant speed below the speed of sound, as time t — )■ oo. 

Our paper is organized as follows: In Section [21 we describe our main result - Theorem 
2.1 - which has four parts. Part (1) is proven in Section 3, parts (2) and (3) in Section 
4, and the proof is concluded in Section 5. Some technical matters are dealt with in two 
appendices. 

Notations: By H^, 1, 2, 3, ■ • • , we denote the Sobolev spaces of complex-valued functions 
on M.^ equipped with the norms 

:= ||(l-A)t/||2. 
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For positive quantities a and 6, the meaning of "a < &" is that there exists a positive constant 
C such that a < Cb. The scalar product of two square-integrable functions, / and g, on 
is given by 

(/, g) ■= j f{x)9{x) dx. 

Acknowledgements 

We are indebted to Daniel Egli, Arick Shao and Israel Michael Sigal for numerous very 
illuminating discussions. Our collaboration on the problems solved in this paper has been 
made possible by a stay at the Institute for Advanced Study in Princeton. We wish to thank 
our colleagues and the staff at the Institute for hospitality. The stay of J.F. at IAS has been 
supported by 'The Fund for Math' and 'The Monell Foundation'. 



2 Statement of the Main Result 

In this section we describe the hypotheses on the potential W and on the choice of initial 
conditions (sec hypotheses {A) and (5), below) and state our main results. 

Theorem 2.1. Suppose the potential W is smooth, spherically symmetric and of rapid decay 
at spatial oo, and suppose that the positive constants M, m and A are of order 1, and the 
constant po > is sufficiently small. 

Then there exists a positive constant eq, such that 

(A) if initially the state of the Bose gas is close to the ground state, (3 = 0, specifically if 
11(1 + |a;p)^t/3o||^3 < eo; and 

(B) if the initial speed of the tracer particle is below the speed of sound in the Bose gas by 
a sufficient margin, specifically if \vo\ = ]g|-Po| is such that \vo\ < \J'^~ ^o, 

the following results hold true: 

(1) For any time t > 0, the speed \vt\ = ^-^ of the tracer particle is strictly below the speed 
of sound, i.e., \vt\ < - C{^o, Po), for some C(eo,po) > 0. 

(2) At large times, the motion of the tracer particle approaches a uniform ballistic (inertial) 

motion: There exists some Poo £ satisfying jj\Poo\ < \J~^> such that Pt — >■ Poo, 
t — >■ oo. 
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(3) In an inertial coordinate system moving with velocity j^Poo, the particle will come to 
rest, i.e., there exists Xoo € M'^ such that Xt — jjPoot — )■ Xoo, as t ^ oo, 

(4) The wave function Pt describing the state of the Bose gas approaches a "traveling wave" 
accompanying the particle, in the sense that there exists a function Poo ^ -L^(M^) with 
the property that 

11(1 + |x - Xt\)-^[Pt - /3oc(- - Xt)]\\oo ^0, ast^oo. 
This theorem is proven in the remainder of our paper. 

As a prehminary to the proof of Theorem 12.11 one must estabhsh local and global well- 
posedness of the equations of motion fll.9p and f ll.lOp . This is easy, because, for an arbitrary 
given particle trajectory {Xt}Q<t<oo^ the equation for [3t, i.e., Eq. f ll.lUp . is linear. Details 
can be found in [2] and will not be repeated here. 

3 An Application of Energy Conservation: Proof of Part (1) 

To prove part (1) of Theorem 2.1, we make use of energy conservation, which holds because 
(11. 9p and (ll.lOp are the equations of motion of an autonomous Hamiltonian system. 

The only non-positive term in the Hamilton functional (I1.12p is l^fp^ /^a W^Ref3 dx. 
This term can be bounded by applying the Schwarz inequality, which yields 

2v^| / RepW^dx\ < \ [ {RePfdx + ^ / {W^fdx 

JrS J^3 a J^3 

= X [ {Rel3fdx + ^ [ W^dx. (3.1) 

Putting this estimate back into the Hamilton functional H (see f ll.l2p ). and using our as- 
sumptions on the initial conditions, we obtain 

H{Xo,Po,l3oJo)<^-^ + 7^ [ \VPo?dx + 2\! \ReP^\^dx+^! W^dx. 
2M 2m V V A 

For any time t, Eqs. fll.l2p and ( 13. ip also yield the lower bound 

H{Xt, Pt, A, A) > !^ - ^ / W'dx. 

These upper and lower bounds, together with energy conservation, i.e., H{Xt, Pt, Pt, Pt) = 
H{Xo, Pq, Po, Po), imply that there is a constant c > such that, for any time t, 

\Pt\' <\Po\' + c{\\Po\\li+Po). 
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Hence if fmax := ig[|-PoP + c(||/3o|li/i + Po)]^ < \/ 2^ then it follows that, for any time t, 



i.e., the speed of the tracer particle remains strictly in the subsonic regime, for all times t. 
This proves Part (1) of Theorem 12. 1[ 

4 Effective Dynamics of the Tracer Particle: Proofs of Parts (2) 



To start with our proofs of parts (2) and (3) of Theorem 12.11 we recast the equations of 
motion in a convenient form. We observe that equation (11.101) is merely real-linear, rather 
than complex-linear, in (3t, and hence it is convenient to rewrite it as a system of equations 
for Re(3t and Im/?^. 

We define a vector function : R'^ — t- by 




(3.2) 



and (3) 



hi(x-Xi) 



lmf3t{x) 



Then and ffLTOD become 



Xt =-Pt, Pt = 







ht) 



(4.1) 



ht =H{t)ht - 





w 



(4.2) 



where H{t) is the 2x2 matrix operator given by 




(4.3) 



To further simplify the equations, we define a new vector function, 6t, by 




(4.4) 
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Rewriting Eqs. (14. ip and (14. 2p in terms of St, we find the equations 







.5t), 
' 

w 



The initial condition for 5 is given by 

5o = -^oH-\Q) 







+ ho, 







H-\t) 







(4.5) 
(4.6) 

(4.7) 
. (To 



and, in the equation for Pt, we have used the fact that 

verify this claim, we determine the explicit expression for H^^{t) and use the fact that W is 
spherically symmetric - the details are straightforward.) 

Applying Duhamel's principle to equation fl4.6p we find that 



5t=?7(t,0)5o + 



M 



U{t,s)H-\s)Ps-V. ds 



--U{t,0)ho~^oU{t,0)H-\0) 



' 


Vp-0 f 


w 


M Jo 





w 

t 



U{t,s)H-\s)Ps-V, ds 





w 



where U{t, s) is the propagator (from time s to time t) generated by the time-dependent 
operator H{-), and we have used the expression for 6o in (14. 7p . 

We now insert Eq. (14. 8 p into the equation of motion f l4.5p for Pt to obtain 



- Po 



M 








f/(t,0)ho) 
,U{t,0)H-\0) 





w 













Uit,s)H-\s)Ps-V, ds 





w 



3, 



(4.9) 



where the terms D^, k = 1,2, 3, correspond to the first, second and third hne, respectively, on 
the right side. These terms depend on the velocity j^Pt, which, by (13. 2p . remains subsonic. 
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In order to have a quantitative measure by how much -o^|-Pf| remains below the speed of 



sound \l -i^-i we introduce the quantity 



M I 

2m' 

\M 

M 



r]:= sup (4.10) 

Y 2m 

The results in Sect. 3 show that < r] < 1. 

The key result of the present section is the following lemma. 

Lemma 4.1. There exists a constant C{ri), which is uniformly bounded on any closed subset 
of [0, 1] not containing 1, such that, 

\D,{t)\<C{ii){l + tr^^,\\{l + x^f\i4n^, (4.11) 

\D2{t)\<C{r^)p,{l + t)-\ (4.12) 

and 

\D:,{t)\<C{r^)p^G{t){l + t)-\ (4.13) 

Here the function G is defined as 

G(t):=max(l + .)3|P,|. (4.14) 

0<s<t 

The proof of this lemma is contained in Appendix |X1 

We now return to equation fl4.9p for P^. By definition, the function G{t) is increasing in 
t, which implies that, for any t > 0, 

G{t) < C{r^)[^4{l + x2)X||h2 + Po + PoG{t)]. (4.15) 

Thus, if Po is so small that C{ri)po < | then 

Git) < 2C{v)[Vp~o\\{1 + x^i^oWn^ + Po] 

which, in view of the definition of G, implies that 

\Pt\<C{l+t)-\ (4.16) 

for any time t > 0, where C is the finite constant given by the right side of the inequality 
above. 
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Next, we show how Parts (2) and (3) of Theorem 2.1 follow from fl4.16p . 

Since, by f l4.16p . Pt is integrable on [0, oo), we conclude that there exists a Poo G such 
that Pt — > Poo, as t — 7- oo. This is Part (2) of Theorem 2.1. 

The position, Yt, of the tracer particle in the inertial coordinate system moving with 
velocity ^Poo is given by 

M 

Its derivative, the velocity of the tracer particle in the moving coordinate system, is obtained 
from 



Yt-^[Pt-P. 



1 

M 



P ds. 



Together with the estimate in f l4.16p . this identity implies that 

\Yt\<C{l + t)-', 

i.e., Yt is integrable in t, which implies that there exists some Y^o such that Yt — )■ Foo, as 
t — 7- oo. This proves Part (3) of Theorem 2.1. 

5 The State of the Bose Gas, as ^ ^ oo: Proof of Part (4) 



In our proof of Part (3) of Theorem 2.1, the function hj 



Re(3t{x + Xt) 
Imf3t{x + Xt) 



has been 





W 



written as a sum of two terms, ^/pqH ^{t) 
establish Part (4) of Theorem 2.1, it suffices to prove that 

||(l + |a;|)-%||oo^O, 



and 6t; see Eq. (14. 4p . Hence, in order to 

(5.1) 



as t -T- oo, with /3oo = y/PoH ^(oo) 





w 



As follows from Eq. (14. 8 p for 6t, we only need to estimate the decay in time of the 
following three vector functions: Ui{t) := f/(t,0)ho, f/2(t) := U{t,0)H-\0) ^ 



W 



and 



U,{t,s):=U{t,s)H-\s) 







In these estimates it is used that the potential W is smooth and of rapid decay at spatial 
infinity. 
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Proposition 5.1. 

11(1 + \x\)-'U,{t)\\^ < (1 + t)-i(l + x^KWh^. (5.2) 
||(l + k|)-if/2(t)||oo<(l + t)-\ (5.3) 

and 

\\(l + \x\)-'Us{t,s)\U<{l + t-s)-\ (5.4) 
This proposition is proven in Appendix B. 

Next, we use Proposition 15 . II to complete our proof of decay estimates on 6t in a weighted 
L°°-space, thus estabhshing Part (4) of Theorem 2.1. Taking the ||(1 + |a;|)~^(-)||oo - norm 
of both sides of Eq. 04.81) and then using Proposition 15. H we conclude that 

2 nt 

11(1 + |a;|)-%l|oo < V 11(1 + \x\)-'Ukmoo + 11(1 + \x\)-'U3{t, s)\UPs\ ds 

k=l "^0 



<il+t)-'+ I (l + t-s)-'\Ps\ ds 







which is the desired estimate. 



A Proof of Lemma 14.11 



In what follows we start by recasting the expressions appearing in Lemma 4.1 in a convenient 
form. Subsequently we sketch the main ideas of the proof. 

A key element of the proof will be to understand the propagator U (t, s) from time s to 
time t generated by H{t). The operator H{t) is the sum of two terms 



H{t) = iJo + § ■ v., 



where Ho is given by 

Ho-.-- 



^A + A) 



2m 
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The observation that these two terms commute enables us to derive an exphcit expression 
for U{t,s): 



U{t,s) = e 



Ho(t-s)(Xt-Xs)-S7^ 



(A.i: 



In our calculations, it is useful to diagonalize the matrix operator Hq and hence the operator 
^Ho{t-s) ^ which, by Eq. (A.l), leads to the diagonahsation of the propagator U{t, s). 



To simplify our notations, we use units such that 

2m = A = 1. 



(A.2) 



We define a 2 x 2 matrix operator, A, by 



A 



zV-A + 1 -V-A + 1 
By Fourier transformation and direct computation we find that 

1 



and 



with L given by 



2i 



^-lg//0(t-s)^ 



^iL{t—s) 







-iL{t-s) 



(A.3) 



Inserting one of the identities A ^A = AA ^ = Id2x2 at appropriate places in the terms 
Di, D2 and of Eq. (14. 9p . we obtain the expressions 

Di{t) = -Re{V^V^W, ie''^e^^'-^'>>^^L'^[tV-A + 1 Re/3o + Im/3o]), (A.4) 



D^it) = -poRe(v^V.W^, e^*^e(^*-^°)-^^[L - ■ d,]-'[-A + 1]-'^W), (A.5) 



and 



(A.6) 
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To prove decay estimates on these terms we use that e'^^^~'^^^e^^*~^^^'^'' and e^iLei^t-Xo)-Va: 
are oscillatory operators in momentum space. More precisely, we Fourier-transform the 
functions inside the scalar products on the right sides of Eqs. (A. 4) - (A. 6), and then employ 
standard integration-by-parts techniques. 

We start with analyzing Di. Without loss of generality, we assume that 

{Xt - Xo) ■ v.. = \Xt - Xo\d,,. (A.7) 
After Fourier transformation, and by introducing polar coordinates, we obtain 

/•27r PIT POO 1^ _^ I 

Di{t)=Re / / e''f^y^+p' — '-i^''"^^p'^Fi{p,e,a) dpdOda 
Jo Jo Jo 

+ Re [ If e''f^'^--^^^"''^^p^F2{p,e,a) dpdOda, 
Jo Jo Jo 

where Fi := sinO W Re/3o niO, a) and Fa := -i^^^nU^M^ n(e, a), and n is the unit vector 

n{6, a) := {cosasinO, sinasinO, cosO)'^ (A. 8) 

Obviously Fi and F2 are complex-valued functions whose smoothness and decay rate at 
p = 00 are determined by those of W, Re/3o and Im/3o. 

To prove an appropriate decay estimate we integrate by parts in the variable p, using the 
recipe 

it + + - ^-r^cosO 

The denominator, a/1 + H — ■f= — — —cos d, on the right side can be seen to be strictly 
and uniformly positive by using the observation that, for 77 < 1 (see (I4.10p ). we have that 

' "Icos < < 1, (A. 10) 



t 

as follows by inserting fl4.10p into = 7I /q § ds\. (We note that the speed of sound, 

2m' ~ because we have assumed that 2m = A = 1.) Consequently 



JlVp' + ^=-\-^i^cos9 1-r/' 
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The factors p^, p^, respectively, appearing in the integrands of the two terms for Di(t) 
enable us to integrate by parts three times, the boundary terms being zero, and we find that 

mt)\ <t-'C mil + \x\'f(3o\\n^, 

where C{r]) comes from controlling the denominator in flA.Op (see also (lA.lip ). It is quite 
tedious, albeit not difficult, to understand how the factor ||(1 + |a;p)^/3o||w2 emerges from 
these calculations. We skip details and refer to a previous paper pj, where a very similar 
problem has been dealt with in a general setting. 

A direct argument shows that 

mt)\<cma+\x\'rMn^, 

which implies that our estimates do not blow up at t = 0. 

The last two estimates clearly imply the desired bound f l4.1ip . 

Bound on D2{t): After Fourier transformation and introduction of polar coordinates, as 
in our analysis of Di{t), we obtain that 

I'2tt p-K poo _^ 

Jo Jo Jo 
The key observation is that the function 

F := [v/rT7 - § ■ fio, a)r\i + p'r'mMd, «), 

with n{6,a) defined in (lA.Sp . is a real-valued, smooth function of fast decay, which follows 
from the fact that, under our hypotheses on the potential W, {Wl"^ is real, spherically sym- 
metric, smooth (in p^) and of fast decay. Thanks to the presence of "Re i" in front of the 
integrations on the right side, we may extend the range of p from [0, oo) to (—00,00) and 
obtain 

1 /-Stt pn poo _^ 

D2{t) = -poRei / / e^'P^^y^+P'-^^^'^'^^p^Fip^e^a) dpdOda. (A.12) 

2 Jo Jo J~oo 

This enables us to use the recipe in flA.9p to integrate by parts as many times as needed, 
which then yields the desired estimate 

\D2{t)\<poCmi + t)-\ (A.13) 
Finally, we turn to the proof of our bound on the term D^, which has the form 

D^{t)=po I K{t,s)-Ps ds. (A. 14) 
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By arguments very similar to those used in the analysis of D2, before f lA.12|) . we find that 

p2-IT PIT POO _^ 

K{t,s) = Re / / / e'(*-Wi+p'-M^--^]p2^2(p',^,«) f^pc^^f^«, 

Jo Jo J -00 



where F2 is a real-valued, smooth function of rapid decay, found after a similar analysis as 
in D2- Here we rotate the axis to make {Xt — X^) ■ 5^ to be \Xt — Xa\dxr^- 

As in our analysis of D2, we integrate by parts three times and then apply a direct 
estimate to obtain 

\K{t,s)\<Cmi + t-s)-\ (A.15) 
Plugging this into (]A.14p we find that 

\D-m <poC{r]) [ {l + t-s)-^\Ps\ ds 
Jo 

<PoCiv) [ (l + t-sr^il + s)--' dsG{t) 
Jo 

with G as in (KW\ . 



B Proof of Proposition 15.11 



The idea underlying our proof is not very difficult. We exploit the dispersive nature of the 
propagators U{t,s) to prove that 6t tends to 0, as t — 00, in a weighted L°°— norm. In 
what follows, we discuss this idea in some detail. 



As in Appendix \^ we assume that 2m = A = 1. 

We start by deriving simple expressions for U{t, s)H~"'{s), n 
defined in flA.SP can be used to diagonalize these three operators: 

f/(t,s)i/-"(s) 



0,1,2. The matrix A 



(B.i: 



--AA-^U{t, s)A[A-^H{s)A]-''A-^ 
--A 









Ki(t,s,n), K2(t,s,n) 
K3(t,s,n), Ki(t,s,n) 



-i{t-s)L {Xt-X,)-d.^ 



A~ 



iB.2) 
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with 



and 



2iV-A + 1 M 

Among all the terms to be analyzed we choose to study the following one 

D ^ ^*^e(^^-^°)-^-it:e^o, (B.3) 

which is a component of the vector function Ui{t) introduced in Sect. 5. All other terms 
can be estimated in a very similar way, the estimates being somewhat easier. 

We propose to prove the estimate 



3 



11(1 + \x\)-'D\\^ < ||(-A + 1)5(1 + \xf)iRePo\\2. (B.4) 
By Fourier transformation we have that 



(27r)'^ Jk3 \k\ 



Introducing polar coordinates we obtain 

f27r /"Tr roo 



D^C [ I I e'Pf^^'''^-''e'^P^^+'^^--^^^^'"^^pF{p,e,a) dpdeda, 
Jo Jo Jo 

for some constant C e C, where F — -\/r+~p^ sin 9 i?e/3o, and f{6, a) e is a real vector. 
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To derive the desired decay estimate we integrate by parts, using ( 1A.9p . and then take 
absolute values. We find that 

1^1 < ^ ' / / / -—[\d,F\{p,e,a) + \F\{p,9,a)]dpd9da 
J Jo Jo Jo + P 



< 



:i + \x\ 



1 + p')-^\dpF\{p,e,a)\U + \\{l + p')-^\F\{p,e,a) 



t 

<^^^^E + l)^^^^e/3o||i + ||(-A + l)ii?e/3o||i] 

k=l 

<^^^ll(-A + l)t(l + |^P)ti?e/3o||2. 

By a direct estimate on D, it is seen that our estimates do not blow up at t = 0: We have 
that 

\D\<\\{-A + l)l{l + \x\')lRef3o\\2. 
Together with the estimate above, this implies flB.4p . which is the desired result. 
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